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A new optical topology and signal readout strategy for a laser interferometer gravitational wave 
detector were proposed recently by Braginsky and Khalili . Their method is based on using a 
nonlinear medium inside a microwave oscillator to detect the gravitational- wave-induced spatial shift 
of the interferometer's standing optical wave. This paper proposes a quantum nondemolition (QND) 
scheme that could be realistically used for such a readout device and discusses a "fundamental" 
sensitivity limit imposed by a higher order optical effect. 
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I. INTRODUCTION AND SUMMARY 



Laser interferometer gravitational wave detectors 
(LIGO, VIRGO, GEO 600, TAMA) are designed to de- 
tect small perturbations h in the spatial metric due to 
gravitational waves (GW) passing through the 
Earth 01 . Being very far from major astrophysical 
sources [H], these detectors are likely to encounter GW's 
that are very weak, so the detectors must be correspond- 
ingly sensitive — e.g., the first LIGO interferometer will 
be able to detect GW's with h ~ 3 x 10~ 21 in the fre- 
quency band of 30 — 300Hz. Improving the sensitivity of 
measurement may be necessary to achieve the first GW 
On ■ detection and will surely be necessary to improve the 
O \ event rate. 

One of the major noise sources in traditional interfer- 
5^ \ ometers is the so-called shot noise. What is being de- 
bJQ' tected is the phase shift of the output optical wave 



WoptTGW^ (1) 



where w Q pt is the angular frequency of the optical wave 



and t g w is the half-period of a gravitational wave. For 
coherent optical pumping the uncertainty in the phase 
due to shot noise is given by A0 = 1 / \J Nqw , where 
A^gw is the number of photons introduced into the inter- 
ferometer during tqw ■ Thus a gravitational wave can be 
detected if 

A^GW > AUn S 77 1 ~r^2- ( 2 ) 

Therefore, in order to increase the gravity-wave sensitiv- 
ity of the interferometer, we have to increase the number 
of photons in the resonator (and hence the consumed 
laser power) as N oc 1/h 2 . On the other hand, the pres- 
ence of the large number of optical photons in the res- 
onator poses severe technical and fundamental 
problems . Among the technical problems are distortion 
of mirrors due to overheating, and large laser power con- 
sumption pi. The fundamental problem is that photons 
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in the interferometer will randomly buffet the mirrors 
inducing random motion indistinguishable from the mo- 
tion produced by a gravitational wave. Balancing this 
radiation pressure noise and the shot noise produces the 
Standard Quantum Limit (SQL) for monitoring the dis- 
placements of the test masses j|. 

Recently Braginsky and Khalili have proposed a new 
way to improve the sensitivity of an interferometric GW 
detector without increasing the interferometer's optical 
power H . Their method entails a new type of GW read- 
out based on a microwave oscillator containing an opti- 
cally nonlinear medium, which is placed inside the GW 
detector's high quality Fabry-Perot resonator. The ad- 
vantage of this readout method is that, unlike conven- 
tional interferometers, it does not require large optical 
power circulating inside the FP resonator in order to 
achieve high sensitivity. In section 2 the principles of 
this scheme are briefly outlined and some numerical es- 
timates are quoted. 

Section 3 describes a potentially practical Quantum 
Nondemolition (QND) strategy which can be used in the 
Braginsky-Khalili readout system (BK-meter) . We show 
in Sec.3 and Appendix A [cf. Eq. (pi])] that a QND mea- 
surement can be performed within a narrow frequency 
band centered around 



= V mcL W ° pt ' 

where w op t and N are the frequency of light and the num- 
ber of photons stored in the FP resonator respectively, 
L is the distance between the end mirrors of the FP res- 
onator and m is the mass of each of the test masses to 
which the mirrors are attached. For N — 2.8 x 10 20 , 
L = 4km, m = 10kg, aj opt = 3 x 10 15 one obtains 
f2o/27r = 60Hz, which is within LIGO band. For the 
resonator's relaxation time of 10 seconds (as assumed in 
||) the necessary laser power to achieve this number of 
photons inside the resonator is ~ 9 Watt. 

We demonstrate in Appendix B that the bandwidth 
Afl of this measurement determines the optimal power 
input to the microwave oscillator : 

Woptimal ~ W-SQL^, (4) 

where PFsql is the power input necessary to achieve the 
Standard Quantum Limit sensitivity at frequency Qq ; cf. 
Eq. ( pi| ) where the expression and the numerical estimate 
for Wsql are given. The signal-to-noise ratio achieved by 
this QND measurement is greater by a factor of i/f^o/ASl 
than the SQL: 



NJ Qm V An Wsql 



(5) 



Section 4 and Appedix B discuss a higher-order optical 
effect in the BK readout system and derive the sensitivity 
limit that it imposes . In particular, thermally excited 
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mechanical modes in the test masses will, after interact- 
ing with light inside the FP resonator, produce a "double 
conversion" of photons, which will be registered as noise 
by the detector; cf. Eq. and Eq. (fffi). 

II. PRINCIPLE OF OPERATION OF THE BK 
METER 

The layout of the BK meter is shown in Fig. 1 (for 
more detail the reader is referred to ||). Three freely sus- 
pended mirrors — A,B and C — form walls of an L-shaped 
Fabry-Perot (FP) resonator which supports a standing 
optical wave, driven by a laser at end A or C. Section A- 
B of the resonator would be in one arm of the LIGO (or 
other) vacuum system, and B-C in the other. The block 
D containing two thin slabs of non-linear medium (Fig. 2) 
is sandwiched between two thin focusing lenses two focal 
lengths apart. The lenses and the block are attached to 
mirror B. 

When the polarization tensor of a gravitational wave is 
aligned with the arms of the FP resonator, the distances 
between A and B and between B and C will change in 
counterphase-i.e. when one is increasing, the other one 
will decrease. This will produce the a net spatial shift 
of the standing optical wave with respect to mirror B, 
thus changing the amplitude of the optical field within 
the two slabs of nonlinear medium. The slabs have cubic 
nonlinearities that are equal in magnitude but opposite 
in sign. They are positioned symmetrically with respect 
to the crest of the standing optical wave as shown in 
Fig. 2. Block D, which contains the slabs, is placed in 
between the plates of a capacitor which in turn is part of 
a microwave oscillator. 

The spatial shift in the optical standing wave produces 
changes of electric field in the first and second slabs that 
are equal in magnitude and opposite in sign: 

8Ei = -SE 2 . (6) 

Since the two slabs have the opposite nonlinearities, 
= > the change in the index of refraction is 

the same for both of them: 

Sn 1 = Sn 2 = 4Trxf ) E 1 5E 1 . (7) 

This change in dielectric constants of the plates in turn 
changes the value of the microwave oscillator's capaci- 
tance, thus producing a shift in its resonant frequency: 

fe c = (8) 

where uj e is the frequency of the microwave oscillator, 
uj opt is the frequency of the optical wave, and 
K = 16n 2 x^ IN huJoptUJe/Vc. Here I is the width of each 
of the nonlinear slabs and V is the volume of the capac- 
itor. This shift is seen as a phase shift in the readout of 
the microwave oscillator: 
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8<f) = Suj e T* = -Kuj op tT*h, 



(9) 



where r* is the oscillator's ringdown time. 

Braginsky and Khalili compare this with the tradi- 
tional optical readout schemes in which the phase shift 
of the recombined optical wave is detected: 

54> ov t = Lo opt T*h, (10) 

where r* is the ringdown time of the two traditional 
FP resonators, one in each arm of the interferometer. 
For x (3) = KT 14 cm 2 /Volt 2 (f US ed silica) , E% pt = 
10 7 Volt 2 /cm 2 (optical breakdown of fused silica) they 
calculate K in Eq. (||) to be of order 1, so for r* pt ~ r* 
the responses of both systems in terms of phaseshift are 
of the same order, 

0opt ~ 4>e- (11) 

For coherent pumping in both cases the uncertainty in 
the phase is A.(f> ~ 1/ VAgw, where A^qw is the number 
of photons (optical in conventional interferometers and 
microwave in the BK readout system) introduced into 
the interferometer during an averaging time (half the GW 
period). So to achieve the same level of sensitivity one 
needs to pump the same number of photons in both cases, 
but the power needed by the BK meter is smaller by 
a factor of w op t/w c ~ 10 4 . The BK estimate for the 
microwave power is 

Wc = ~ IWatt (12) 

for N c ~ 10 20 . For more detailed estimates the reader is 
refered to Ref. @ . 



III. QND FOR THE BK READOUT SYSTEM 

Any readout system that monitors the displacement 
of the mirrors must exert on them a fluctuating back 
action force, thus enforcing the Heisenberg uncertainty 
relation . As a consequence of this, all straightforward 
displacement measurements run into the Standard Quan- 
tum Limit (SQL) [|, § 

Ax SQL = J—, (13) 

V TTL 

where Ax is the minimal uncertainty in displacement of 
a free mass m monitored over a time interval r. This 
SQL for displacement can also be written in terms of 
the limiting spectral density of the mirrors' displacement 
fluctuations j|: 

^ SQL (^) = 4)2 ( 14 ) 

mi r 

where Q is the frequency. Then Ax = VS^ QL Ail, and 
for fl ~ AS1 ~ 1/t one recovers Eq. (O) 
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The SQL for a free mass is by no means a fundamental 
limit; it can be overcome by a variety of techniques || 
which are known collectively as Quantum Nondemolition 
(QND) measurements. All previously proposed QND 
schemes that are applicable for conventional GW inter- 
ferometers utilize highly non-classical states of light, and 
none of them are practical because of technical difficul- 
ties (most especially because of the large required optical 
pumping power and because losses so easily destroy the 
non-classical states of light). In this paper a different 
strategy is proposed, one which does not require the de- 
liberate creation or detection of any non-classical state 
of light and thus can be more practically implemented. 
This scheme, however, is confined to narrow-band mea- 
surements. 

We begin by describing the backaction mechanism by 
which the BK readout system enforces the Heisenberg 
Uncertainty relation on the measurement of the test-mass 
position. The quantum state of the BK microwave oscil- 
lator satisfies the usual phase-number uncertainty rela- 
tion 

Acb c AN c > ~. (15) 

The more accurately the BK- meter reads out e , the 
larger will be the fluctuations AiV e in the oscillator's 
number of microwave photons. The nonlinearity will 
transform AN C into an uncertainty of the optical index 
of refraction of the slabs: 

16tt 2 y (3) ?L 

S ni = -5n 2 = £ SN C , (16) 

eV 

where SN e is the fluctuation in N e , Sni and Sn 2 are the 
resulting fluctuations in n\ and n 2 , e is the coefficient of 
dielectric permitivity and V is the volume of the capaci- 
tor. Braginsky and Khalili have argued Q that Sni and 
Sn 2 cause a redistribution of the optical energy between 
the left and the right parts of the FP resonator, thereby 
giving rise to a net difference in the forces buffeting the 
mirrors: 

SF = K <^^ (17) 

where L is the total length of the FP resonator. This 
fluctuating force will cause fluctuations in the positions 
of the mirrors, thus causing fluctuations in the spatial 
shift of the optical field with respect to the mirror B and 
the nonlinear slabs attached to it: 



where <5i(0), SF(fi) and SN e (fl) are Fourier components 
of the corresponding quantities, m is the mass of each of 
the mirrors , and the factor | comes about when motion 
of all three mirrors is taken into account. 

Now we are ready to describe our QND method, but 
first the following simple remark must be made. Suppose 
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for a moment that all of the mirrors are rigidly fixed. As 
already mentioned above, fluctuations in N c , by changing 
the optical coefficient of refraction of the slabs [Eq. (Eq)], 
will redistribute optical energy between the left and right 
parts of the FP resonator. A straightforward calculation 



[Eq. ( A16 ) in Appendix A] shows that this alone will 



change the optical field inside the nonlinear slabs so that 

SEi = -6E 2 = Eo^^Sn, (19) 
cV2 

thereby simulating a spatial shift of the optical field as in 
Eq. (H). Here I is the width of the nonlinear slab , n = 
n\ = n 2 , 5n = dni = —5n 2 and E is the peak amplitude 
of the optical standing wave inside the FP resonator. 

Now if we release the mirrors, the back action ( jlq ) will 
affect our reading as well, and the total fluctuations of 
the optical field inside the nonlinear slabs will be given 
by 

SEx(Q) - SE 2 (n) = 




where r = L/c; cf. Eq. (A17) of Appendix A. From the 



above equation we see that for a given frequency O = S^ 
we can adjust N in such a way that SEi (SI) = 5E 2 (CI) = 
and thus the readout system does not register any fluc- 
tuations due to the back action (but only for that value 
of SI). Thus a QND measurement is performed. The re- 
lationship between the QND angular frequency Slo and 
the number N of optical photons in the Fabry-Perot res- 
onator is 



see Eq. (g) of Sec. I. 

The essential reason that this readout is QND is that 
it registers not only the fluctuations of the mirrors' dis- 
placement x due to backaction [the second term in large 
parentheses in Eq. (17)], but also directly the back-action 
force (the first term). Thus a position-momentum cor- 
relation is introduced into the measurement procedure, 
and such correlations are known to make QND possible 
{§. For L = 4km, Q/2tt = 60Hz, w opt = 3 x lO^s" 1 , 
and m = 10kg the necessary number of optical photons 
to perform QND is 

N=- Q } m f T ~ 2.8 x 10 20 . (22) 

The QND measurement described above is clearly nar- 
rowband . But one can dynamically tune the frequency 
at which the QND is performed by changing the laser 
power and thus changing the number of optical photons 
N in the resonator, provided that the frequency of the 
signal changes slowly compared to the ring-down rate of 
the optical resonator. 



G 



Appendix B considers a particular scheme for mea- 
suring of the phase of the microwave oscillator. In 
this scheme the oscillator is coupled to a transmission 
line, and the physically measured quantity is the phase 
quadrature of the outgoing electromagnetic wave propa- 
gating along the transmission line. Having specified fully 
the measurement model, we find that if the bandwidth 
of measurement is An then the signal-to-noise ratio for 
the narrow-band QND measurement can be as high as 

^ ) QND W SQL ^ ^ ' 

i n r n "+ An m n 2 L 2 \h(n)\ 2 

= —— — - / ail. 

27rAllJ no _ AQ h 

The above signal-to-noise ratio is achieved when the 
pumping power of the microwave oscillator is given by 



^optimal 2 



v 2 ( n V f L \ 1 n o 
I j hujr An 



32w 2 x {3) Nt V^opt 

where V is the volume of the capacitor and Wsql is the 
minimal power necessary to achieve the SQL sensitivity 



level; cf. Eq. (|B13|) of Appendix B and Eq. <M of the 
introduction. For V = (.01mm) 3 , uj e — lO 1 ^ -1 and 
for other parameters having numerical values as in Eq. 
(||), we get W^ptimai = O.lWatt (n /An) (cf. the lkW 
of the optical power required to achieve the SQL in a 
conventional interferometric scheme). 

While we have not devised a general proof, it seems 
likely that no other microwave readout scheme can oper- 
ate with a power less than in Eq. (|24|); expression J24|) 
is probably a general relation for optimally disigned mi- 
crowave readout schemes. 



IV. HIGHER-ORDER OPTICAL EFFECTS; 
FUNDAMENTAL SENSITIVITY LIMIT 

In this section we identify and discuss a fundamental 
limit on the sensitivity of the BK readout system — a limit 
that applies whether or not the system is being operated 
in a QND mode. 

In an interferometric GW detector, mirrors are in- 
stalled on the surfaces of test masses, which have inter- 
nal elastic mechanical modes of frequencies ^l m j1n > 
12 kHz. The noise curve of the interferometer will have 
large peaks near these frequencies. When photons of fre- 
quency w pt interact with walls oscillating with the fre- 
quency n m , some of the photons will be up or down 
converted to frequencies w pt ± n m . These up or down 
converted photons in turn interact with the "noisy" walls, 
and if there is a non-zero component of the mirrors' mo- 
tion at n m ±n, then some of the photons will up or down 
convert a second time to frequencies wo ± n. If n is the 
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frequency of detection, then this second-order process of 
double frequency conversion will be registered by the BK 
readout system as a signal from a gravitational wave. 

The perturbation theory for FP resonators with mov- 
ing walls is worked out in detail in App endi x C. Here 
just the main result is quoted. From Eq. ( |C26| ) the noise 
curve in units of 1/ -s/Hz is given by 

1 ^opt&f 



where T e is the temperature of the test masses, j m is 
the damping rate of their mechanical modes, and e ~ 

max(5x (3) /x (3) >^A)- Here 5 X = X? + X ( i\ X (3) = 
Xi i ^ is t ne wavelength of light in the resonator and 
SI is the spatial offset of the slabs from the position in 
shown Fig. 2. For L = 4 km, fi = 60rad/sec, T e = 300K, 
ft m = 7.2 x 10 4 rad/sec, m — 10kg and 7,„ = 10 -8 fl m we 
get the noise level of 



~ 10 -28 V(fir + e)/VHz (26) 

It is not unimaginable that future interferometers will 
achieve sensitivities y/Eh ~ 10 _29 /\/Hz for low frequency 
(10 — 100Hz) narrow band signals (by e.g. using the QND 
technique described in this paper). In this case, Eqs. J25|) 
and ( p6| ) show that higher order effects will give rise to a 
"fundamental" low frequency noise limit of magnitude 

V^^O-^^pVHz. (27) 



V. CONCLUSIONS 



In this paper we have shown that a practical QND 
measurement might be possible for a narrow-band mea- 
surement by a gravitational wave interferometer using a 
BK readout system. Also it was shown that second-order 
effects set a "fundumental" limit on the precision of the 
measurement. 



VI. ACKNOWLEDGEMENTS 



I want to thank Vladimir Braginsky for suggesting the 
problem, for interesting and useful discussions, and for 
wonderful Pasadena ice-cream, Kip Thorne for reading 
carefully over the manuscript and making many useful 
suggestions, and Farid Khalili and Sergei Vyatchanin 
for interesting discussions and excellent tea. I also 
thank Glenn Sobermann for his advice on the prose, 
and Michele Pickerell and Ruben Krasnopol'skiy for their 
help with computer graphics. This work was supported 
in part by NSF grant PHY-9424337. 



8 



APPENDIX A: PHYSICS OF THE NONLINEAR 
MEDIUM INSIDE THE BK FABRY-PEROT 
RESONATOR 



First consider one slab of nonlinear medium positioned 
inside a FP resonator. Let xi (Fig. 3) be the total path 
length from the left mirror to the left edge of the slab, 
X3 be the path length from the right mirror to the right 
edge of the slab and I be the width of the slab. For sim- 
plicity of the calculation we assume I <C A where A is the 
wavelength of light in the resonator. Also for convenience 
define t\ — Xi/c, r 2 = l/c, t 3 = x 3 /c. 

The eigenfrequencies u> of this optical resonator were 
worked out in ||. They satisfy the following eigenequa- 
tion: 



sin(wr) — (n — 1) sin(nu;T2) 

(Al) 



sin (ijJTi) sin (u;r 3 ) H — cos (ustA cos (ustA 
n 



where r = T\ + r 3 + 11T2 ■ This equation has approximate 
solutions 

us = us + {n 2 - l)u) ^- {cos [us (n - t 3 )] - 1} (A2) 

where ujq = irk/r, and k is any integer. When the slab's 
index of refraction n changes, us changes accordingly: 

dus nui°T 2 , r o .1 -> 

— = — — {cos [us (n - r 3 )J - 1} (A3) 

The total optical energy contained in the resonator is 

U = Nhus (A4) 

where N is the number of optical photons . We can find 
all of the forces acting on the mirrors by taking deriva- 
tives of U with respect to t\ and t 3 . For example, 

Nh dus , A x 

F left = — A5 

C OTi 

and 

Nh dus , . 

bright = 7^ A6 

c <9t 3 

where Fi e ft and F r ; g ht are the forces acting on the left and 
the right mirrors respectively, with the positive direction 
being out of the resonator. When taking derivatives of us 
one has to keep in mind that usq also depends on t\ and 

The force acting on the slab of nonlinear medium is 
-Fieft — -Fright • The total spatial shift of the optical wave 
with respect to the slab due to the forces acting on the 
end mirrors and the slab itself is 



5X(Q) = "2^2 [ F ^tm - bright (fi)J , (A7) 
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where "tildas" stand for Fourier Transforms. If i<i e ft and 
-Flight are produced by a fluctuating inde x o f refrac tion 
n = n + Sn then on substituting Eqs. ( |A5| ) and ( |A6| ) 
into Eq. (A7) we get 



3Nh ( d 



d \ duj 



2mfl 2 \ dri 9t3 J dn 



Sn{n). 



(A8) 



By then putting Eqs. (A3) and (A8) together we obtain 
3Nhn uj 2 t 2 



n 2 



■sin[w(n -T 3 )]Sn(Cl). (A9) 



For the two slabs of opposite nonlinearities (Sni — 
—Sn 2 = Sn) in the configuration of Fig. 2 their two con- 
tributions add up to give 



mc iV t 



(A10) 



The above expression is a manifestation of the back ac- 
tion as explained in Sec. II. 

Now let the amplitude of the optical electric field in 
the left part of the resonator be 



E, 



left 



E sin — 

V c / 



(All) 



where x is the spatial coordinate with the origin at the 
left wall. Then the field in the middle of the left slab is 
given by 



Ei — Eq 
+0 



1 (nUT2 

sin (ujti) H — cos (uiTi) sin I — - — 
n V 2 



Now 



dE 1 dE dui 

— — ~ — — Sin(WTi) + EqTi cos(wri) — . 

dn dn dn 



(A12) 



(A13) 



In the case when two slabs are present inside the FP 
resonator, their contributions to the frequency and field 
changes add up linearly (since the perturbations are very 
small) . For the configuration of Fig. 2 we see from 
Eq. (|A3|) that du/dn — 0, so 



dE x 
dn 



dEo 
dn 



sin(wri). 



(A14) 



But F le ft OC Eq, so 

dEo I dn 1 AFicft / dn 



-Flcft 



t d du> 
uj dn dr\ 



-3T7^- (A15) 



Putting Eq.(A3) and Eq.(Al4) into Eq.(A15) and doing 
exactly the same calculation for the second slab, we ob- 
tain 



5Ei = -5E 2 = E 



->opt 



y/2c 



Sn 



(A16) 
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for the case when mirror B is in the middle of the 
reso nato r. Combining this with the back action from 
Eq.([A10|) we finally get Eq.(§0|) of Sec. 3: 



5E! (O) - 6E 2 (Q) 



V2nE ^ (1-6 




opt 



APPENDIX B: CALCULATION OF OPTIMAL 
MICROWAVE POWER AND SIGNAL-TO-NOISE 
RATIO FOR A QND MEASUREMENT 

Consider the microwave oscillator as shown on Fig. 2. 
In order to get information about the phase of the oscilla- 
tor, we have to couple it to the outside world. Whatever 
the nature of this coupling is, it will cause dissipation 
of the induced oscillations and hence, by the fluctuation- 
dissipation theorem, give birth to a fluctuating compo- 
nent of the oscillator's current. 

For concreteness, we model this coupling by an open 
transmission line of impedance R. We assume that the 
oscillator, consisting of the capacitor C and inductor £, 
is driven on resonance by a generator G with a voltage 
output of amplitude Vo (see Fig. 2). We also assume 
the transmission line encompasses all of the dissipation 
present in the oscillator, i.e., more generally, that we can 
access all of the information escaping from the oscilla- 
tor. And finally, we set the temperature of the outside 
world to (in reality, one will have to cool the oscilla- 
tor to temperatures below the ones corresponding to a 
microwave frequency) . The ingoing vacuum modes drive 
fluctuations in the circuit as described above, and the 
phase of the outgoing wave contains information about 
the phase of the oscillator. 

The ingoing modes are described by the positive fre- 
quency part of a voltage operator 

V in = / duV ' Rhua in {Lu)e~ luJ \ (Bl) 
Jo 

where a\ n {uj) is the annihilation operator for the 
ingoing mode of frequency to normalized so that 
(0|ai n (o;)a i ' n ([j')|0) = 5(uj — ui'). Then the fourier compo- 
nent of the outgoing wave is 

-a + ifi 



V out (u e + fi) = \J RTiuj e 7^a in (uj e + Q) 

a — i\l 

+ (B2) 
2\l + la 

where a — R/C is the ringdown rate of the microwave 
oscillator and 5ui e is the variation in the oscillator's res- 
onant frequency due to fluctuating optical fields in the 
slabs of nonlinear medium, as explained in Sec. II: 
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8uj r 



8n 2 x i3) lNhu> opt uJe 8Ek - SE 2 



V2VL 



En 



(B3) 



Here V is the volume of the capacitor. The change of the 
optical field inside the slabs is given by 



8E x {n) - SE 2 {n) = V2E nojoptT 2 [1-6 



TLC 2 T f2 2 



5n(Q) 



+V2E ^x s (fl) 



(B4) 



where the first term on the right hand side is due to 
the fluctuating index of refraction of the slabs [cf. Eq. 
( pp| ) of Sec. Ill and discussion therein], and the second 
term is due to the GW-induced relative displacement x s 
of the slabs with respect to the standing optical wave. 
The fluctuations Sn of the indices of refraction of the 
nonlinear slabs in the above expression are caused by the 
voltage fluctuations on the plates of the capacitor, which 
in turn can be traced to the incoming vacuum modes of 
the transmission line: 



Sn(n) 



2Tr X (3) VMuT e V uj 2 e 
1 a{2tt + ia)d 2 



ain(^e + fi) + a\ n (uj e - tt) 
(B5) 



where d is the distance between the plates of the capac- 
itor. Collecting Equations @), @, @ and @ 
together, we can write down the expression for the phase 
quadrature of the outgoing wave in the transmission line, 
which is the measured readout signal: 



r(n) = [VoutK + n) - 

y/Rhuj e a + id 



Vo a — itt 
S^Vox^INHlo 2 ^ 



a in (u t 
optL ° {x t 



- fi) 
Q) - 



/V 



n) 



(B6) 



V (2Q + ia)VLc 
8tt 2 x {3) VoVM^I 



1 



6NH < pt 



RV(2n + la) 
The measured x is then given by 

8T7 2 X (3) V^eWl 



n 2 



a in (u> e +Sl) + aL (w e - £1) 



^measured (^) — ^s(^) 



V(2Q + tot) 
VLc(2n + ia)y/huT e 



1 - 



6Nh w 2 opt 



16w 2 X ( 3) lNhuj 2 pt uj e VW 



mc 2 T fl 2 
a in (u e + O) - a\ n (w e - fi) 



am (w e + 0) + a\ n (ui e - fi) 



(B7) 



where W — V 2 /R is the power pumped into the mi- 
crowave oscillator by the generator G. The correspond- 
ing spectral density of the Gaussian noise seen by the 
readout system is 



s x (Q) 



8ir 2 x {3) l 



V 



huj e w 



AVt 2 



VLc 



1 6NH <g 
mc 2 T n 2 

(An 2 + a 2 )huj e 

w 



+ 



(B8) 
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The first term on the right hand side corresponds to the 
back-action noise and the second term corresponds to the 
intrinsic noise of the measuring device. 

We aim to perform a measurement with a narrow fre- 
quency band centered around the frequency fto at which 
the back-action noise is zero: 



/6iV^ (BQ) 



We write S x as a Taylor expansion in frequency around 
ft : 

\2 



_(ft-_ft ) 2 , B(n ) 
i 2 

'n 



s x (si)~A(n )w " ft2 u/ 1 (B10) 



where A and B can be read from Eq. (B8). If the relevant 
bandwidth is Aft then we place a limit (ft — fto) 2 < Aft 2 
and 

Aft 2 B(Qn) 
S x (Sl) < A(ft )W-^- + -L^. (Bll) 

Minimizing the right hand side of the above equation 
with respect to W , we find the expression for the mini- 
mum noise in a fixed bandwidth: 

Sxoptimali^) < 2VAB — — = (B12) 
1 A 2 rAft Aft 

-2^1v""ft7~^ SQL(i2oJ ft^' 

which is achieved at the input power 

V 2 ( ft \ 2 /L\ 2 1 ft 



Woptunal 3 27r 2 x (3) 2 7V T ^opt J \, I J Jw e Aft ' 

(B13) 

In the above expressions A is the wavelength of light 
inside the FP resonator and SsQL(fto) is the Standard 
Quantum Limit noise at the frequency ft for a free mass. 
Clearly, the signal-to-noise ratio for this narrow-band 
measurement is ^/fto/Aft greater than that in the case 
of the SQL : 



S_ 

N 



APPENDIX C: PERTURBATION THEORY FOR 
FABRY-PEROT CAVITY WITH MOVING WALLS 

In this appendix we derive a formal series for the opti- 
cal field inside a Fabry-Perot resonator which is pumped 
by a monochromatic laser beam and the walls of which 
are free to perform motions small compared to the wave- 
length of light A. The expansion parameter is 5x/X, 
where Sx is the change of length of the resonator. For 
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our purposes we are only interested in expanding up to 
(Sx/ A) 2 ; and we use this formal series to derive Eq. (25) 



The following situation is considered: for simplicity 
we assume that light is pumped on resonance by a laser 
beam E m = ae~ 

M"optt-kx) though trie left mirror which 
is at rest and has reflectivity r and transmissivity T . For 
concreteness it is assumed that the fluctuations in length 
Sx originate from the motion of the right mirror which 
is assumed to be perfectly reflecting. Further, we as- 
sume that the plain wave approximation is applicable 
and hence the optical field inside the resonator satisfies 
the onc-dimcnsional wave equation: 

c 2 ^-)A(x 1 t)=0 (CI) 



dt 2 dx 2 , 
The general solution of of the above equation is 

A(x,t)=.f(t+-)+g(t--) (C2) 



where / and g are arbitrary functions. The boundary 
conditions at the left mirror (x = 0) and at the right 
mirror (x = L + <5x) read respectively 

g(t) - rf(t) = Tae-^ (C3) 

and 

f(t + T O + ^j+g(t-r - S -^)=0 (C4) 

where To = L/c. Eliminating g(t) from these two equa- 
tions, we get 

/(* + to + ^) -rf(t-n>-?f\ ( C5 ) 

= Tae~ luJ ° ptt e l ^° ptT °e l ~^ Sx 
or, expanding in Sx up to second order, 

f(t + r ) - rf(t - t ) = Tae- tu °** t e w °* n (l + i^Sx - ^-Sx 2 ^ 

--[f(t + T )+rf(t-r )]Sx (C6) 
c 

We take a Fourier tranform of the above equation and 
then solve it by iterations: 

/(«) = + fM(w) + fW(u) + ... (C7) 

where 

f^(uj) = l£Ls{uj-u opt), (C8) 
1 — r 

x 2z^ opt ra _ 

7 M " c[(l-r)cos(a;T )-z(l + r)sin(a;To)](l-r) dn ' J W ° ptj ' 

(C9) 
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/«(u,) = ^ tTa 



c 2 [(1 — r) cos (wro) — i{l + r) sin (ojtq)] (1 — r) 

, cos{u)'t )5x(uj' - ui op t)Sx(uj - uj') /r , inN 

ClLO 7- r -. ; > 7- r — ; -. ; r. (OIUI 

(1 — r) cos(w'tq) — 1(1 + r) sin (a/ To) 

When writing down the above terms we took into account 
the fact that 1 — r <C 1. The structure of is clear: 
it corresponds to upconversion of light at frequency uj pt 
to an intermediate frequency u>' and then from lu' t o u , 
with u>' being integrated over. From Eqs. ( |C2| ) and ( |C^ ) 

A(x,u) ~ -2«sin^af) /(a;). (Cll) 

The BK readout system detects the square of the am- 
plitude of the optical field: 

S(x,t) = \A(x,t)\ 2 (C12) 
= S (0) (x, t) + S^\x, t) + S( 2 \x, t) + ... 

where 

SW(x, Q) = 4sin 2 (^z) C 2 <5(f>), (C13) 



(C14) 



and 



s (2) {x,n) = - 




( 2w op tC 
V c 



oj pt+0 — SI 



+ S i n (2zf2. x } 



(C15) 



£(n)£(n') 



(fo(fi')(fo(fi - tt')dtt'. 



In the above expression C = Ta/(1 — r) and 
£(Q) = (1 - r) cos(r2r) + (1 + r) sin(fh"). 

In real interferometers Sx represents, for example, mo- 
tion of the surface of the mirror due to the thermal exci- 
tation of the test mass's internal modes. In what follows 
the contribution from the internal mode of lowest fre- 
quency is considered and then it will be shown that the 
sum of contributions of all the higher modes will have the 
same order of magnitude. It is assumed that the ther- 
mal noise is a Markoff Gaussian process, and therefore is 
described by the following equation: 



F(0) 



n 2 -n 2 



(C16) 



where VL m is the eigenfrequency of the mechanical mode, 
7 m is the damping rate and F(TL) is the Langevin force 
satisfying 



(F(n x )F(n 2 )) 



D 
2^ 



<5(fi 1 + fi 2 ). 



(C17) 



where D = kBT e ~/ m /m* is the velocity diffusion rate. 
Here m* is the effective mass of the mode (approximately 
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given by the mirror mass to), ks is Boltzmann's constant 
and T e is the temperature of the enviroment. To calcu- 
late the spectral density of the fluctuations of (the 
goal of this analysis) we will need the 4-point correlation 
function : 

D 2 

(F(Sl 1 )F(Sl 2 )F(Sl 3 )F(Sl 4 )) = ^-2 ^("i + ^)S(Sl 3 + Sl A ) 

07T 



(C18) 



+s(n 1 + si 3 )S(si 2 + st 4 ) + d(si 1 + si 4 )6(si 2 + si 3 )} 



Using the above expression, Eq. (|C15| ) and Eq. ( |C16| ) 
we obtain 

(5 (2) (x, fii)5 (2) (x, fia)) = {2uj Q C f [MiS (Six) 5 (fi 2 ) + Ah (x, Sli) S (fii + Sl 2 )] . 

(C19) 

Here M 2 (Sl), which characterizes the spectral density of 
fluctuations of 5^ 2 \ is given by 



D 2 f 

M 2 (x, Sl) = — [K (SI, SI') K (-SI, —SI ) + (C20) 
K (SI, SI') K (-SI, -SI + SI')} dSl', 



where 
K(Sl,Sl') 



C(Sl')(Sl' 2 - Sl 2 n + i lm SV) \ (Sl - SI') 2 - Sl 2 n + i lm (SI - SI') 



sm 



£(fi-Q') 



(C21) 



sin (^x) 



£(fi) 



It is possible to integrate Eq. ( |C2C| ) exactly, but it is 
clear that the main contribution will come from mechan- 
ical and optical resonances, SI' — Sl opt and SI' = Sl rn . 
For 7 m <C (1 — r)/r (which is the c ase f or, e.g., fused 
silica) the major contribution in Eq. ( C2C| ) is due to the 
mechanical resonances: 

D 2 

M 2 (x, SI) ~ — (2u ovt C) 2 [Kt (x, SI) + K 2 (x, SI) + K 3 (x, SI)} 



(C22) 



where 



K 4 (x,Sl) 



2sin 2 ( "opt-n^ gin 2 ^o Pt +n m ^ + sin 4 ^o P t+n m ^ + gin 4 ^o Pt -n m ^ 

16sin 4 (n m r)^fi 2 7 m 
(C23) 



K 2 (x,Sl) 



16sin i (fi m T)(fir) 2 rj^^ 2 7„ 



-, (C24) 



l(i 



, in 2 ^opt+n m ^ +gin 2 ^ opt -n m ^j ^2(5^^ 



16sin' 3 (fi m T)(r!T)r!4 i 02 7m 

(C25) 



We are only interested in detection frequencies such that 
CIt <C 1. Then for the configuration of nonlinear slabs 
shown in Fig. 2 the main contribution to the noise in 
the BK meter readout will come from K% and K3. The 
spectral density of the displacement noise will be 

1 / k B T e \ 

^Opt "TV? 



s '-^{^) ^a5#w (fe + '>' (C26) 

where e characterizes the degree of positioning error and 
the mismatch of nonlinearities of the two slabs: 



(C27) 



Here 61 is the spatial offset of the central point between 
the two slabs. 




Equation ( |C2q ) is the main result of this Appendix. 
It's implications are discussed at the end of Sec. IV. 
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See Fig. 2 




Fig. 1 

The BK readout system. 

Two thin lenses focuse light on a block D containing 
nonlinear medium. Block D is shown on Fig. 2. 




Fig. 2. 

This is an enlarged view of the block D from Fig. 1. 

Two slabs of nonlinear medium are positioned at the points 
of maximal gradient of the intensity of the optical standing wave. 
The spacial shift of the optical standing wave changes the resonant 
frequency of the microwave LC oscillator. 




Fig. 3 

The changing coefficient of refraction of theslab of nonlinear medium 
will redistribute optical energy between the left and the right parts 
of the Fabry-Perot resonator. 



